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Abstract: We extend the design of radially symmetric three-dimensional 
invisibility cloaks through transformation optics JT] to cloaks with a surface 
of revolution. We derive the expression of the transformation matrix and 
show that one of its eigenvalues vanishes on the inner boundary of the 
cloaks, while the other two remain strictly positive and bounded. The 
validity of our approach is confirmed by finite edge-elements computations 
for a non-convex cloak of varying thickness. 
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1. Introduction 



It was recently found that transformation optics open new avenues in electromagnetic cloaking, 
either through their heterogeneous anisotropic effective material parameters (transformation 
optics, (HlHl) or through low index materials U or negative refractive index materials (4). In- 
terestingly, the invisibility is preserved in the case of an intense near field Q, when the ray op- 
tics picture breaks down. The mathematics behind the scene have been known from researchers 
working in the area of inverse conductivity problems J6). The first experimental realization of 
an invisibility cloak, chiefly achieved in the microwave regime JT), suggests that cloaking will 
be limited to a very narrow range of frequencies. However, it will not be perfect since the cloak 
is necessarily dissipative and dispersive, and some of its tensor components are singular on the 
inner boundary. The latter drawback can be overcome by considering a generalized transform 
lU [9] in an upper dimensional space and then projecting the resulting metric on the physical 
space, and this leads to non-singular tensors of permittivity and permeability. Alternatively, one 
can design an approximate structured cloak via homogenization lfT0lfTTlfT2l . although the rapid 
growth of the field, fueled by a keen interest of the optics community, promises a large panel of 
new technological applications. 

A non trivial question to ask is whether one can design cloaks of non-spherical shapes. A 
parameterization of the cloak's boundaries was proposed by three of us to design cylindrical 
cloaks of an arbitrary cross-section ifTTl . Its extension to the general three-dimensional case 
now requires to parameterize the inner and outer boundaries of the cloak as some surfaces with 
varying radii 

p(0,0)=ao.o+ {fl m ,«cos(m0 +n(j>) +b mj , sm(m9 +n<j>)} , (1) 

(m,n)GN 2 \{(0,0)} 

which is somewhat out of reach running the COMSOL multiphysics package on a 256 Gb RAM 
computer: Compared to a spherical cloak, the tetrahedral mesh needs be further refined due to 
the complexity of surfaces involved. However, to demonstrate the versatility of our Fourier- 
based approach, it is enough to analyse arbitrary cloaks built with surfaces of revolution i.e. 
generated by rotating a curve about an axis. 

In the present paper, we discuss the design of axially invariant three-dimensional cloaks 
with an arbitrary cross-section described by two functions R\(<p) and /?2(0) giving an angle 
dependent distance from the origin. These functions correspond respectively to the interior 
and exterior boundary of the cloak. We shall only assume that these two boundaries can be 
represented by a differentiable function. Their finite Fourier expansions are thought in the form 
RjiQ) = a o + Tin=i a Q n cos ( n 0)' J = l'^' where p can be a small integer, and a J Qn = for 
n > p for computational easiness. Note that our approach encompasses the case of toroidal 
cloaks |[T3l but cloaks with irregular boundaries such as cubes, fall beyond the scope of this 
study. Nevertheless, contrarily to previous three dimensional numerical studies lfl3l [121 . our 
scheme can be used in the design of cloaks with non-convex boundaries. 

To illustrate our methodology, we compute the electromagnetic field diffracted by a cloak 
with rotational symmetry about the z-axis. We perform full-wave finite element simulations in 
the commercial package COMSOL, when the cloak is illuminated by an approximate plane 
wave (generated by a constant electric field on the upper surface of the computational domain). 



2. Change of coordinates and pullbacks from optical to physical space 



2.1. Material properties of the heterogeneous anisotropic cloak 

The geometric transformation which maps the field within the full domain p < ^2(0) onto the 
annular domain R\ (0) < p' < RiiQ) can be expressed as: 

e' = e , < < 27r , (2) 

</>' = , -n/2 < <j> < k/2 

where < p < RiiQ)- Note that the transformation maps the field for p > R^iQ) onto itself 
through the identity transformation. 

This change of co-ordinates is characterized by the transformation of the differentials through 
the Jacobian: 

3(p(p',f),M) 

This change of coordinates amounts to replacing a homogeneous isotropic medium with 
scalar permittivity and permeability e and fi, by a metamaterial described by anisotropic het- 
erogeneous matrices of permittivity and permeability given by ll5l fT4ll 

j^=£T-\ and = /iT -1 , (4) 

where T = J r J/ det(J) is a representation of the metric tensor in the so called stretched radial 
coordinates. Note that there is no change in the impedance of the media since the permittivity 
and permeability undergo the same transformation. 
After some elementary algebra, we find that 
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where 



and 



Cl3( * )= * 2( *Wf)-«i(f)) 2 w + * l(f) wW)F rff ' () 

for Ri((j>') < p' <RiW)- Elsewhere, T _1 reduces to the identity matrix (en = 1, C13 =0 and 
p = p' fov p' > R 2 {(j)')). 

2.2. Singularity analysis of the transformation matrix 

To exemplify the symmetric nature of the coefficients of T~ 1 , and due to the role played by its 
first entry in the following singularity analysis, we show the variation of (T ) u in Fig.Q] We 
note that this coefficient varies between zero and three. The fine mesh on the inner boundary 
of the cloak is also apparent from Fig.Q] as are the non-convex and symmetric features of the 
cloak. 



In the case of a spherical cloak, C13 vanishes and T 1 reduces to Diag( t p ^ ,cn,cii). We 
note that the first eigenvalue vanishes on the inner boundary and the other two remain constant, 
which is consistent with the singularity analysis led in |[T5l , This is unlike the circular cylindri- 
cal case whereby one eigenvalue goes to zero while the other one goes to infinity on the inner 
boundary ||5l, 





Fig. 1. 3D plot of (T ')[[, as given by Eq. [5] Eq. [6] and Eq. [7] within the cloak with 
boundaries given by Eq.[lO] The symmetries are noted. 

However, in our case, the cloak is of an arbitrary shape, and it is therefore illuminating to 
look at the behaviour of its permittivity and permeability tensors's eigenvalues. The eigenvalues 
of (0 are found to be 



1 - ^+p 2 +jiP a , (-iy./f '?3+p a +jiP fl v 
a j- w2 +—y{ ) 
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;' = 1,2, and A 3 = cn 



(8) 



We can therefore see that Xj, j = 1,2 are spatially varying functions of p'(p, 0), such that Aj = 
when p = i.e. at the inner boundary, while A2 > 0. Importantly, A3 is a strictly positive constant 
for a given angle 0' (i.e. a function independent upon p). Last, we checked that when there is no 
longer a symmetry of revolution about one axis, all three eigenvalues are also spatially varying 
with 0', but A3 remains independent upon p. This reflects the fact that the geodesies for light 
within the arbitrary cloak follow more and more complex trajectories when we perturb the 
geometry away from the spherical cloak. 

3. Finite elements computations 

We would now like to further investigate the electromagnetic response of the cloak to an inci- 
dent plane wave from above, see Fig. [2] For this, we choose the electric field E as the unknown: 
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-k 2 e'E = 



(9) 



where k = co^jioEo = co/c is the wavenumber, c being the speed of light in vacuum, and s!_ 
and \i' are defined by Eqs. ©. Also, E = E, ■ + E c /, where E, is the incident field and E^ is 
the diffracted field which satisfies the usual outgoing wave conditions (to ensure existence and 
uniqueness of the solution). The weak formulation associated with Eq. |9] is discretised using 
second order finite edge elements (or Whitney forms) which behave nicely under geometric 
transforms (pull-back properties) 03). 
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of wavenumber k = 2^/0.3 incident from above on an non-convex invisibility cloak. 



For the sake of illustration, cf. Fig. [2] let us consider a cloak with inner and outer boundaries 
expressed as 

Ri(<j)) =O.2 + O.O2cos(40) , R 2 (<j>) = 0.4 + O.O2cos(80) . (10) 
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Fig. 3. 2D plot of *AEj +i?2 +£3 generated by a slice of Fig.[2]in the xz-plane for y = 0. 



We report the computations for the magnitude of the total electric field in Fig. [2] (3D plot), 
Fig. [3] (2D plot in the xz-plane) and Fig. H] (2D plot in the yz-plane) for a plane wave incident 
from above at wavenumber k = 27T/0.3 (units are in inverse of a length, say fim for nearly 
visible light (UV)). 

Around 4.10 5 elements were used in this computation, which corresponds to about 2.7 10 6 
degrees of freedom. While the convergence of the numerical scheme has been checked by 
considering different types of meshes, the large size of the system means we were not able to 
further refine the mesh for the computational resources at hand. 
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4. Conclusion 

In conclusion, we have proposed a design of an arbitrarily shaped cloak using a Fourier ap- 
proach. We only assumed that the cloak displays a symmetry of revolution about one axis in 
order to reduce the computational complexity of the problem (extension to completely arbitrar- 
ily shaped cloaks is a straightforward matter). Cloaking has been confirmed numerically for 
an incident plane wave in resonance with the concealed region and an analysis of the cloak's 
singularity has been carried out. 
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